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ABSTRACT 

■ Decomposing the shear signal into E and B-modes properly, i.e. without leakage of B-modes 

\ into the E-mode signal and vice versa, has been a long-standing problem in weak gravita- 

tional lensing. At the two-point level this problem was resolved by developing the so-called 
ring statistics, and later the COSEBIs; however, extending these concepts to the three-point 
level is far from trivial. Currently used methods to decompose three-point shear correlation 
functions (3PCFs) into E- and B-modes require knowledge of the 3PCF down to arbitrary 
£j . small scales. This implies that the 3PCF needs to be modeled on scales smaller than the min- 

imum separation of 2 galaxies and subsequently will be biased towards the model, or, in the 
absence of a model, the statistics is affected by E/B-mode leakage (or mixing). 
In this paper we derive a new third-order E/B-mode statistic that performs the decomposi- 
tion using the 3PCF only on a finite interval, and thereby is free of any E/B-mode leakage 
| while at the same time relying solely on information from the data. In addition, we relate 

I/"") . this third-order ring statistics to the convergence field, thereby enabling a fast and convenient 

calculation of this statistic from numerical simulations. We note that our new statistics should 
be applicable to corresponding E/B-mode separation problems in the CMB polarization field. 
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j> ! 1 INTRODUCTION 
• i-H ■ 

■ Cosmic shear, the distortion of light from distant galaxies by the tidal gravitational field of the intervening large-scale structure, is an ex- 
$H , cellent tool to probe the matter distribution in the universe. The statistics of the image distortions are related to the statistical properties of 
5^ ■ the large-scale matter distribution and the geometry of the universe, and can thereby be used to constrain cosmology. Current results already 

demonstrate the power of cosmic shear observations at constraining the clustering amplitude cr 8 and the matter density Q m (e.g., Fu et al. 
2008; Schrabback et al. 2010; Huff et al. 2011). Furthermore, cosmic shear provides an ideal tool to study dark energy through measuring 
the growth of structure with large future surveys like KIDS 1 , DES 2 , LSST 3 (Ivezic et al. 2008), or Euclid 4 (Laureijs et al. 2011). The large 
volume probed by these surveys will enable us to measure not only the power spectrum, but also higher-order statistics with unprecedented 
precision. As the evolved density field is non-Gaussian, the three-point correlation function and its Fourier space equivalent, the bispectrum, 
contain significant cosmological information complementary to the more commonly used two-point statistics and are a powerful tool for 
breaking parameter degeneracies (Takada & Jain 2004). 

The upcoming weak lensing experiments will limit the statistical uncertainties to the percent level. In order to extract cosmological informa- 
tion from these cosmic shear experiments, the increased data quality needs to be accompanied by a thorough treatment of a wide range of 
systematic errors, from photometric redshifts and galaxy shape measurements to the removal of astrophysical contaminants. 
If the shear estimated from observed galaxy shapes is solely caused by gravitational lensing, then it should consist only of a "gradient com- 
ponent", the so-called E-mode shear. B-modes (or curl components) cannot be generated by gravitational light deflection to leading order, 
and higher-order corrections are expected to be very small. Hence observing any B-mode pattern indicates remaining systematics in the shear 
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analysis. 

Decomposing the observed shear field directly into E/B-modes (e.g Bunn et al. 2003) is complicated by the complex mask geometry of weak 
lensing observations. At the two-point statistics level, an E/B-mode decomposition is commonly performed using the aperture mass disper- 
sion (Schneider et al. 1998) and related measures (e.g. Crittenden et al. 2002), which can be calculated from the measured shear two-point 
correlation function (2PCF) and is thus not affected by the masking geometry. However, these methods assume that the 2PCF is known either 
from 6 = to some finite angular value (aperture mass dispersion) or to arbitrarily large separations. However, in reality the 2PCF can only be 
measured on a finite interval [# min , 6 m . llL ], where the lower boundary is caused by inability to measure the shape of image pairs with very small 
angular separation. As Kilbinger et al. (2006) pointed out, lack of shear-correlation measurements on small scales leads to an underestimation 
of the aperture mass dispersion on small scales and causes an apparent mixing of E- and B-modes with this type of estimator. Schneider & 
Kilbinger (2007); Eifler et al. (2010) and Schneider et al. (2010) develop statistical measures for an exact E/B-mode decomposition based on 
2PCFs known only on a finite interval [9 min , # max ]- 

At the three-point statistics level, Jarvis et al. (2004) and Schneider et al. (2005) introduced E/B-mode separating shear measures which as- 
sume knowledge of the 3PCF down to arbitrarily small scales. Shi et al. (201 1) derived a general condition for the E/B-mode decomposition 
of lensing three-point statistics, but the construction of filter functions with finite support based on this condition is far from straight forward. 
In this paper we derive an extension of the 2PCF ring statistics (Schneider & Kilbinger 2007; Eifler et al. 2010) to an exact E/B-mode de- 
composition of shear three-point correlation functions on a finite interval. 

In order to constrain cosmology with third-order shear statistics it is important to obtain the corresponding predictions from a large suite of 
cosmological numerical simulations in a reasonable time and with limited computational effort. We facilitate this by giving an expression of 
the third-order ring statistics in terms of the convergence field, thereby avoiding the time-consuming calculation of the shear 3PCF for each 
simulation. 



2 SHEAR THREE-POINT CORRELATION FUNCTIONS 

We first introduce the shear three-point correlation function (3PCF): Consider a triangle in the complex plane with vertices X; and let y^X;), 
H = 1, 2 be the Cartesian components of the shear at point X,. Unless otherwise noted, we will assume that the triangle is oriented such that 
Xi, X 2 , X 3 are ordered counterclockwise around the triangle. We define x t = Xj - X 3 and x 2 = X 2 - X 3 to be the sides of this triangle (c.f. 
Fig. 1). We will use x, to refer to complex numbers or vectors interchangeably, and denote their magnitude as x,. 
The Cartesian components of the shear 3PCF are defined as 

Wxi.x 2 ) = (y^(x 1 )r,(x 2 ) r/t (x 3 )} , (i) 

where we have assumed that the shear field is statistically homogeneous so that y^x depends only on the side vectors x,-. Since one cannot 
form a tri-linear scalar from the product of three shears, the behavior of the Cartesian components of the shear 3PCF under rotations is 
complicated. In order to write the 3PCF in terms of tangential (y t ) and cross components (y x ) of the shear which are parity eigenstates 
and have relatively simple transformation properties, one can project the complex Cartesian shear y c = y l + iy 2 into tangential and cross 
component with respect to a chosen direction a, with polar angle a h 

y(X,; a,) = y t (X,; a,) + iy x (X,; a,) = - [y,(X,) + iy 2 (X,)] e" 2 ^ = -/(X^" 2 ^ = -/(X,) a' 2 /a 2 . (2) 

If the directions of projection a, are defined in terms of the vertices X; and thus do not depend on an external coordinate system, then the 
tangential and cross shear are invariant under rotations of the triangle (Schneider & Lombardi 2003; Takada & Jain 2003; Zaldarriaga & 
Scoccimarro 2003), and the 3PCF of these shear projections will only depend on the side lengths x, and the orientation of the triangle (clock- 
wise or counterclockwise). In the following we will use the centroid projection, where the shear at vertex X; is projected along the direction 
q, connecting X, with the centroid X = (X t + X 2 + X 3 )/3, and a t is the polar angle of this projection direction (see Fig. 1 for an illustration). 
Following Schneider & Lombardi (2003) we define the (complex) natural components of the 3PCF which have relatively simple transforma- 
tion properties 

r<°W 2 ,* 3 ) ee <y(X 1 ;a,)y(X 2 ;a 2 )y(X 3 ;a 3 )> = - <y c (X,)y c (X 2 )y c (X 3 )>e- 2l <">™> 
I< 1 \xuX2,x 3 ) = <y*(X 1 ;a 1 )y(X 2 ;a 2 )y(X 3 ;a 3 )> = - <y c * (X,)/ (X 2 )y c (X 3 )> e^'^" 3 ' 

I< 2 \ Xl ,x 2 ,x 3 ) = <y(X 1 ;a 1 )y*(X 2 ;a 2 )y(X 3 ;a 3 )>, r (3) (x 1 ,x 2 ,x 3 ) = <y(X 1 ;a 1 )y(X 2 ;a 2 )y*(X 3 ;a 3 )> . (3) 

r (0) is invariant under cyclic permutations of arguments; the other three components transform into each other: r (1) (xi , x 2 , x 3 ) = r (2) (x 3 , x\ , x 2 ) = 
r (3) (x 2 , x 3 , Xi), etc.. A different parameterization of oriented triangles is in terms of two sides and their inner angle, e.g. xi,x 2 , and <f> (c.f. 
Fig. 1). We choose the convention <p e such that <p > corresponds to Xj, X 2 , X 3 being ordered counter clock wise ("positive 

orientation") and <p < corresponds to clock wise ordering ("negative orientation"). 
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3 E/B-MODE SEPARATION 

To construct integrals which separate third-order E- and B-mode correlations we start from the circle statistics C{9) (Crittenden et al. 2002; 
Schneider et al. 2002) which geometrically separates E- and B-modes by measuring the mean tangential and cross component of the shear 
on a circle of radius 9 around the origin 

i r 2n i r 2n 

C{6) = C<(0) + iC x (0) = — J o d^ (y, + iy x ) (9, if,) = -— dip e" 2i y (9, f ) , (4) 

where tf, is the polar angle on the circle, and in the last step we have rotated the tangential and radial shear into the cartesian components. 
Following Schneider & Kilbinger (2007) we now consider the shear inside an annulus &\ < 6 ^ & 2 and define the ring statistics R 

r'h 

R = R t + iR x = \ dBW(fi;&i,& 2 )C(.0), (5) 
which is a function of two radii and &j, and where W(8; i?,, &j) is a normalized weight function 



60 W(9;&i,&j) = 1, (6) 

and W = outside the annulus, i.e. if < or > From this definition we construct the third-order ring statistics as the correlation of 
the weighted mean shear in three concentric annuli with radii i?i < 6>i, < # 2 < #3 < 9 2 < #4 < &s < O3 < #6 (cf. Fig. 1), 

(RRR)(&) = f 2 de^^&u&i) f ' d9 2 W(e 2 ;& 3 ,& 4 ) f 6 d9 3 W(9r,& 5 ,& 6 ) (Cie^^)^)} (7) 

J#] Jl? 3 J# 5 

dfliW(fli;#i,tf4) d<W(0 2 ;# 3 ,#4) de 3 W(e 3 ;& 5 ,& 6 ) (C (0i)C(0 2 )C(0 3 )> (8) 

!>! Jlf 3 Jlf 5 

where we have used = i? 6 ) to denote a six-tuple of radii. Expanding these correlators in terms of the mean tangential and cross shear 

yields 

(KRK) = - (R x R x R t ) - (R x R t R x ) - (R t R x R x ) + i (- (R X R X R X ) + (R x %R t ) + (W) + WWU)] W (9) 

<^W5> (*) = [(%%%) + (R X R X %) + <Kx*tKx> - <K,K x Kx> + i ((R X R X R X ) - (R x %%) + (%R X %) + <*W*x»] (*) (10) 

Note that the imaginary parts of (9, 10) vanish in the absence of parity-violating modes. 

We analogously define the correlators (RR'R) and {R.RR*} and separate E- and B- modes via 

= ^Re [{RRR} + (R'RR) + (RR'R) + <«RR*>] (1?) (11) 
(ft E «£)(#) = ^Re [-3 {<RWR) + (R?<RR) + (WR"R) + (RRJ?y\(fl) (12) 
(ft£« B )(0) = |Im [3 (ftftR) + {R*R3i) + <KR*K> + {RWR*)] (&) (13) 

= ilm [- <K^> + (K'Kft) + <ftR*K> + <7W?K*>] (0), (14) 

where Eq. (11) corresponds to pure E-mode correlations, and Eq. (14) to parity violating third-order B-mode correlations. Equation (12) is a 
parity invariant correlation between E- and B-modes, and Eq. (13) is a parity violating correlation between E- and B-modes. 
For brevity, the mixed terms (12,13) are generalized expressions which are sensitive to B-modes in any of the annuli, i.e (li E 'R^j('ff) = 
({R t R x R x ) + (R x R t R x ) + {R x R x Rt))/3, etc.. Instead one can also consider more localized B-mode measures like 

{R x R x R t )(&) = ^Re [{-RRR) + {R*RR) + {RR*R)-{RRK )}(&), (15) 
which picks up correlations with B modes in the innermost and middle annulus, but is insensitive to B-modes in the outer annulus. 



4 THIRD-ORDER RING STATISTICS 

In this section we derive computationally advantageous expressions for the third-order ring statistics in terms of the shear 3PCF, and show 
their relation to the convergence bispectrum. 

4.1 Relation to the shear three-point functions 

We rewrite the third-order ring statistics in terms of the shear 3PCF by starting from the definition Eq. (7) 
{RRR)(&) = - f 2 d0iW(0i;0i,02) f * d(9 2 t¥(<? 2 ; i? 3 , i? 4 ) f * de 3 W(9 3 ;& 5 ,& 6 ) 

J#l Jl> 3 Jl?5 

x r it r r ^ <?° ^ ^ ■ ^ 
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Figure 1. Left: Concept of the third-order ring statistics (UKR) (d\ , d 2 , $6)- E/B-modes are separated on a finite interval by correlating the shear of galaxy 
triplets located within three concentric annuli, as illustrated by the thick black triangle. The minimum separation between galaxies in the above geometry is 
min(i?3 — d 2 , #5 — #4) and the maximum separation is + #4. 

Right: Geometry of a triangle in the third-order ring statistics. The triangle vertices Xj are located on three concentric circles of radius 9j and have polar angles 
tjjj. The triangle centroid is X. In the centroid projection the shear at each Xy is projected onto the centroid along direction q^, the line connecting Xj with the 
centroid. ips is the orientation angle of vector q^. Finally, is the inner angle of the triangle at X3 which we will use when parameterizing a triangle in terms 
of two side lengths x\ , x 2 and angle <p. 



Noting that X ; = 8j exp(ifa) and using Eq. (3), this can be rewritten as 

x exp(2i(a 1 + a 2 + a 3 -fa-fa- fa)) r m (Xi - X 3 , X 2 - X 3 ) , 



(17) 



where r <0) is the shear 3PCF measured relative to the centroid, so that the a, are the directions of the point X, to the centroid X = (X[ + 
X2 + X3)/3. Owing to circular symmetry, we can set if/ 3 = 0; equivalently, one can use relative polar angles Afa = if/j - if/ 3 and show that the 
integrand depends only on these relative angles. 

As T is measured within discrete angular bins, while the weight functions and geometric factors in Eq. (17) can be evaluated continuously, 
it is numerically more stable to rewrite the third-order ring statistics such that only the three outermost integrals contain the shear 3PCF and 
the inner integrals can be evaluated numerically to arbitrary precision. With = X ; - - X3 = Bj^i - 8 3 for j = 1,2, 

(KRR)(0) = -31 J i2xi J d2 -*2 r< °W x 2) J d03 W(er,&5,&6)W(\x 1 + 3 \;& 1 ,& 2 )W(}x2 + 6 3 \;& 3 ,&4) 
1 

x l 7T7\ 7T* ex P(2i(«i + a 2 + a 3 - fa ~ fa) , (18) 

|Xi + 3 ||x 2 + 3 \ 

where we have used 6 3 to denote a complex number with zero imaginary part for consistency. We have e 2 ""j = q ; /q*, with qi = (2xi - X2)/3, 
q 2 = (2x 2 - xi)/3, q 3 = -(xi + x 2 )/3, and €*) = (x.j + 3 )/6j, so that 



(19) 



Thus, 
CRKK) (0) 



(2. 

X 



— J d 2 Xl J d 2 x 2 r <0) ( Xl ,x 2 ) J d8 3 W(8 3 ;& 5 ,& 6 )W(\x l +0 3 \;& l ,& 2 )W(.\x 2 + 3 \;& 3 ,& 4 ) 

qiqiq 3 x l +o 3 x\ + 0\ 



1 



(20) 



|xj + 3 ||x 2 + 6 3 \ q;q 2 qj x, + 3 x 2 + 3 
Finally, if ^>,- is the polar angle of x,, and <f> = <p 2 - ip\ is the angle between x 2 and Xi, we obtain 

<KKR>(0) = j djcixj j dx 2 x 2 j d0r <o) (x,,x2,0) 

/C -r ^ 1 aia9Q3 xt + 6 3 xt + B 3 

d9 3 W(8 3 ;& 5 ,& 6 ) d Vl W(\ Xl e m +8 3 \;& l ,& 2 )W(\x 2 e' (m+ ' l > ) +8 3 \;& 3 ,& 4 ) - — — ^H-^— 
J |xi + 3 ||x 2 + 03 1 q;q 2 q3 Xi + 9 3 x 2 + 8 3 

= t^tj f dxiXi J dx 2 x 2 J dc/,T {0) (x i ,x 2 ,(f>)Z Q (x 1 ,x 2 , l f>,&), (21) 
where we have defined the complex filter function Zo of the ring statistics in the last step. Note that the ratio of the q's does not depend on 8 3 
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Figure 2. Contours of the filter functions Zo,i as a function of angular scales x\ e [#5 - #2. #6 + #2] and jc 2 6 [#5 — $4, $6 + #4] for = (1', 2', 3', 4', 5', 6') 
with (f> = 71/8 (top row) or <f> = n/4 (bottom row). Contour lines are evenly spaced with separation 0.5 dex ranging from |Z| = 10~ 5 to |Z| = 10 -1 ' 5 , dashed lines 
indicate regions where Z is negative. The filter functions vanish if the triangle configuration (x, ,X2,(f) is not allowed in the ring statistics geometry (c.f. Fig 1). 

and thus the evaluation of the filter function can be further simplified by reversing the-order of integration and moving this factor to the outer 
(if,-) integral. 

Expressions for the other correlations required for E/B-mode separation, which contain a complex conjugate ring statistic H*, are derived 
analogously. For the correlation involving the complex conjugate shear at vertex X ; , the resulting expression analogous to Eq. (21) contains 
instead of r <0) , q ; /q* is replaced by its complex conjugate (corresponding to a, — > —aj in Eq. (3)), and for j = 1,2 the factor x^/x* is 
also replaced by its complex conjugate (corresponding to tf/j — > -if/j in the equivalent of Eq. (16)), e.g. 

CR"KR){§) = J dx lXl J dx 2 x 2 J d0r <1) (x 1 ,x 2 ,0)Z 1 (x 1 ,x 2 ,0,^), (22) 

with 

Ziixuxi,*,*) = f M 3 W(9 3 ;& 5 ,& 6 ) f dip, W(\x,^ + 9 3 \;&,,& 2 ) W(\x 2 ^ + ® + 9 3 \;& 3 ,& 4 ) -J — q ' q ^ *i + ft ^±^ (23) 

J J |xi + 3 ||x 2 + 9 3 \ qiq 2 q; xj + 9 3 x 2 + 6 3 

In the computation of the ring statistics one can choose any (normalized) radial weight function W that fulfills W(0; & 2 ) = even if &, = 
(as the separation in tangential/cross shear is ill-defined on circle of radius 6 = 0). To be specific, we choose 

W(6; &j) = 30(6 - &i) 2 ($j - dfK&j - & t f , (24) 

as in the computation of the second-order ring statistics (Schneider & Kilbinger 2007). The shape of the third-order ring statistics filter 
functions Z ,i based on this choice for W is illustrated in Fig. 2. 

4.2 Relation to the bispectrum 

In order to rewrite the third-order ring statistics in terms of the bispectrum we first relate it to the lensing convergence field k, which is easier 
to express in terms of the convergence bispectrum than the shear 3PCF (see Schneider et al. 2005,for details) as it contains fewer oscillatory 
phase factors. Expressing the ring statistics in terms of the convergence field also speeds up the measurement of in simulations consid- 
erably, as described below. 
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Figure 3. Filter functions associated with K for different ring radii ($,, The left panel shows the radial filter function of the ring statistics, W(9;&i, #;); 
the corresponding compensated aperture mass filter function Ufi h § j {8) is shown in the middle panel, and the right panel illustrates the corresponding Fourier 
transform Ug^^t). 



Consider the convergence field smoothed with a radially symmetric filter U#(ff) with characteristic scale &. If U#{&) is a compensated filter 
f dd 6 U§(9) = 0, this convolution can be expressed in terms of the shear field as 

J d 2 6' UMff\) K(ff) = J d 2 Q 9 (\ff\) y t (0') , (25) 

where U and Q are related by (Kaiser et al. 1995; Schneider 1996; Schneider & Kilbinger 2007) 

2 C 6 f °° 2dff 

Qm = ¥ 1 MV'M-UoW and u *^ = J g —Qd9')-Qd0). (26) 

As shown in Schneider & Kilbinger (2007) the definition of the ring statistics H (Eq. (5)) is equivalent to an aperture mass M^(& h &j) with 
two characteristic scales if 

^w = -W^- (27) 

As the relation between the filter Q and U does not depend on the shape of Q, we can calculate the corresponding compensated filter U# it #.(0) 
as in Eq. (26). The left and middle panel of Fig. 3 show the ring statistics filter W(9; & h &j) and the corresponding aperture mass filter 
for different choices of ring radii As expected from Eq. (26), U is constant for 9 < & h then becomes negative, and is zero for 9 > &j. 

Based on Eqs. (25,26,27), the third-order ring statistics of a pure E-mode field can be computed directly from simulated convergence maps 
by convolving the convergence field with different filters U$ l j )j and correlating three filtered maps. With this approach one does not need to 
calculate the shear 3PCF, which are computationally expensive (e.g. Jarvis et al. 2004). 

Expressing K E as the convolution of k and U^ u § } also enables us to write down he third-order ring statistics of a pure E-mode field in terms 
of the convergence bispectrum B K (l ly l 2 , h) (c.f. Schneider et al. 2005) 

(ftl) W = J dh /, J dl 2 1 2 J d(f> B K (luh, + -2/ 1 / 2 cos0)f7,, 1 ,,, 2 (/,)(/,,3,,, 4 (/ 2 )(/,, 5 ,,, (5 ( A //2 + l 2_ 2 / 1 / 2 cos0) , (28) 

with the Fourier transformed filter function 0(1) = fd99 J (W)U(6). The bispectrum filter functions for the third-order ring statistics are 
illustrated in the right panel of Fig. 3. 



5 CONCLUSION 

Upcoming lensing surveys will provide data of unprecedented quality and enable us to conduct robust measurements of cosmic shear beyond 
the two-point level. These higher-order statistics contribute substantial information to cosmological constraints by breaking parameter degen- 
eracies when combined with second-order shear statistics. Furthermore, three-pt statistics have the potential to improve our understanding 
of systematics effects in the data, e.g. a detection of third-order B -modes can be an additional indicator for unsolved problems in the data 
analysis. 

When extracting third-order information from a high-quality data set it is therefore essential to use robust and unbiased theoretical methods 
that meet the quality of the data. We have introduced the third-order ring statistics, which separates the shear 3PCF into third-order E/B- 
mode correlations on a finite interval [6( min , max ]. Hence this statistic does not require knowledge of the 3PCF down to zero lag, where it is 
impossible to measure. Thus, unlike the third-order aperture mass statistics, it is not affected by apparent E/B-mode mixing (Kilbinger et al. 
2006). 

Our main results are Eqs. (21, 22, 23), which give compact expressions for the third-order ring statistics in terms of the shear 3PCF. Fur- 
thermore, in Sect. 4.2 we give convenient expressions for computing the E-mode ring statistics from numerical simulations, and from the 
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convergence bispectrum which facilitate the comparison with theoretical models for weak lensing three-point statistics (e.g. Valageas et al. 
2011). 

In addition to the cosmological information contained in the E-mode signal, our expression for third-order B-mode correlations opens a new 
window to detect remaining systematics in the data. For example, the various permutations of CReRb.'R.e) allow for an association of B-modes 
with a specific angular scales. 

For the analysis of future shear 3PCF measurements, we recommend using Eq. (1 1) to obtain a clean third-order E-mode signal, and Eqs. (12, 
15) to test for remaining B-mode correlations. 
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